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We discuss systems of ions in segmented linear Paul traps for the quantum simulation of spin models which 
exhibit long-distance entanglement (LDE) in the ground state. The axial potential can be tailored to create 
different coupling patterns suitable to create LDE. Microwave pulses are used to implement various kinds of 
spin-spin interaction Hamiltonians, which can be varied adiabatically thereby realizing the adiabatic preparation 
of the ground states of quantum spin Hamiltonians exhibiting LDE. 
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I. INTRODUCTION 



Entanglement is a central resource for quantum technologi- 
cal applications lj],|2[]. Great effort has been devoted to the 
generation and distribution of entanglement over long dis- 
tances, with the intent of realizing large-scale protocols of 
quantum information and communication technologies [J- 



Particularly intriguing in this context is the prediction that 
certain spin models are naturally endowed with peculiar en- 
tanglement properties in their ground state which could be 
profitable for quantum communication purposes, for example, 
between different spatial regions within a quantum proces- 
sor. Specifically, the concept of long-distance entanglement 
(LDE) has been introduced and discussed in order to identify 
the occurrence of sizeable nonlocal quantum correlations be- 
tween distant, non-directly interacting spins in quantum spin 
chains and networks Il7l - l22ll . This phenomenon emerges in 
models with non degenerate ground states, when the end spins 
(spins at the boundary of the system) interact weakly with 
their immediate neighbors, such that a strongly correlated bulk 
mediates effective interactions between the distant, non di- 
rectly interacting, end spins. In this work we discuss the fea- 
sibility of schemes for the experimental observation of this 
effect using trapped ions as quantum simulators of quantum 
spin models. 

Trapped ions are highly versatile systems which have been 
proven to be extremely effective in quantum technological ap- 
plications. The simulation of quantum models of strongly in- 
teracting quantum matter with trapped ions holds promise for 
the investigation of those quantum dynamics that remain so far 
unexplored due their inescapable complexity [23]. Indeed, the 
natural many-body dynamics of trapped atoms is very rich and 
interesting by itself; on the other hand, in the present work we 
will be mainly concerned with the subtle and intriguing task of 
realizing models that are not directly provided by the natural, 
i.e. non engineered, physics of trapped ions. Although spin 
interactions emerge quite naturally in ion chain systems, engi- 
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neering and control of a desired complex Hamiltonian can be 
highly challenging tasks with high pay-off. S pec tacular proof 
of principle experimental demonstrations II24h27|| have shown 
the potential of trapped ion based quantum simulators. How- 
ever, so far none of these experiments has explored the ground 
state of spin models which are expected to exhibit highly non- 
classical properties. Moreover, no experiment so far has faced 
the simultaneous occurrence of all the following features com- 
bined together: 1) sophisticated shaping of the trap potentials 
in order to suppress the effect of long-range interactions; 2) 
highly controlled adiabatic processes driving the system to the 
ground state; 3) full implementation of Trotterization (Trotter 
expansion) in order to generate the relevant spin-spin interac- 
tions in all the needed directions and components. 

In the following we will discuss schemes to simulate 
ground states of complex quantum systems, exhibiting highly 
nonclassical properties, by the use of trapped ions featuring 
aspects l)-3) listed above. Specifically, we will explore the 
capabilities of trapped ion systems for the quantum simula- 
tion of spin models, and we will apply them to propose the 
first demonstration of LDE in quantum spin chains. LDE is 
a global nonclassical effect which, on the other hand, can be 
monitored by the analysis of only two spins, namely the end 
spins of the chain. It is therefore a sufficiently simple, yet rich 
phenomenon which is ideal to be demonstrated using an ion 
trap quantum simulator. Differing from the previous experi- 
ments cited above J25l - t27ll in which the spin coherent manipu- 
lation is realized with laser fields, here we focus on segmented 
ion traps in the presence of a magnetic gradient where the en- 
gineering of the quantum dynamics is realized by microwave 
fields. 

The paper is organized as follows. In Sec. [TT] we introduce 
the systems and discuss the basic features of the scheme that 
we implement for the simulation of long distance entangle- 
ment. In Sec. [HI] we discuss how to tailor the spin-spin inter- 
actions and we describe the scheme of pulses for the simula- 
tion of spin Hamiltonians. In Sec. [TV] we discuss the results 
for the adiabatic preparation of the ground state and discuss 
the experimental feasibility of the protocol. Finally, in Sec.[Vl 
we draw conclusions and discuss possible outlooks. 



II. THE SYSTEM 

Doppler cooled ions held in a segmented ion trap 
and exposed to a magnetic field gradient realize effective spin- 
1/2 models H30l — [3311 . The effective spin-spin interactions in- 
duced by the magnetic field are of Ising type and can be ad- 
justed by tailoring the axial trapping potential. In particular, 
if the ions are sufficiently cold, such that the ion motion can 
be neglected (the validity of this approximation is discussed 
in Sec. lIVBl l, the effective system of N spins is described by 
the Ising Hamiltonian 
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where the resonance frequencies of the atomic spins a>j de- 
pend on the external magnetic field B(xoj) at the equilibrium 
position of the ion xqj IT34ll . The spin-spin couplings are in 
general long range and their magnitude depends on the trap- 
ping potential and on the spatial derivative of the spin reso- 
nance frequency that, in turn, is determined by the magnetic 
field gradient. They are given by 
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Here, A, whose elements are A^ = 



d 2 V(xi,-x N )\ 



, is the 

x N ) 
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Hessian matrix of the potential energy function V{x\, ■ 
that confine the ions with Xj indicating the position of the ion 
j. Besides, the magnetic gradient allows also for the address- 
ing of individual spins with microwave field which can, there- 
fore, be used to manipulate the spin dynamics. 



A. General considerations 

Spin Hamiltonians with non-trivial ground state correla- 
tions (as in the case of LDE) are in general characterized by 
non-commuting spin-spin interaction terms. This is not the 
case of the simple Ising Hamiltonian (Q]) in which only in- 
teractions terms of the form crcr^ are present. Therefore, the 
simulation of LDE requires the ability to engineer interactions 
along a different axes, described for example by a term of the 
form crfcrf . Such an effective interaction can be induced using 

J k b 

a sequence of tt/2 microwave pules that realize the transfor- 
mation &~ mcr i' cr.e 1Kcr j' - cr x over all the spins. In particular, 
a free evolution sandwiched by two trains of n/2 pulses (each 
pulse addresses a particular ion j = 1,...,/V) with opposite 
phase performs the following transformation 
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and realizes an Ising interaction along the x-axes. In order for 
this transformation to be effective, the duration of the pulses 



have to be sufficiently short so that the evolution due to the 
spin-spin interactions can be neglected. This is achieved with 
a sufficiently strong microwave driving field resulting in a 
Rabi frequency Q » Jj^, On the other hand each microwave 
pulse should operate on a single spin, and its effect on the 
other spins should be negligible. This imposes a limit on the 
maximum allowed intensity of the driving field O <k A, where 
A indicates the frequency difference between neighboring spin 
resonances. 

The simultaneous interaction along z and x can be simulated 
by Trotterization, namely by repeated, fast application of the 
two kind of interactions 1 35]. Provided that the interaction 
time TJn is sufficiently small, it is possible to approximate 
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and to generate a stroboscopic evolution which simulates a 
Hamiltonian that is the sum of two Ising Hamiltonians with 
interactions along two orthogonal axes. 

We also note that typically the parameters in the Hamilto- 



Ising 



nian H'Z;„„ defined in Eq. (JTJ are such that the spin-spin cou- 



pling strengths are much smaller then the single site energy, 
Jjjt <K a>c, which hence dominate the dynamics of this model. 
Nevertheless we note that we are interested in the situation in 
which the system is driven by a series of microwave pulses. 
In this case, as demonstrated in the next section, the relevant 
dynamics is that obtained in a reference frame rotating at the 
driving field frequency. In this representation, the relevant sin- 
gle site energy is in fact given by the detuning h = 0)j- Vj be- 
tween spin resonance frequency and driving field frequency, 
which can therefore be adjusted and controlled during the dy- 
namics. 

These results can eventually be used for the adiabatic prepa- 
ration of the ground state of, for example, XX Hamiltoni- 
ans. The system is prepared initially in the ground state 
of a sufficiently simple Hamiltonian which is easy to pre- 
pare: In our case it consists of the ferromagnetic/fully polar- 
ized spin state which is the ground state of the Ising Hamil- 
tonian with finite magnetic field (i/;„,v/ fl / = fi h/2 £ ; - cr. - 
h/2 Tjj,k Jj.k^j ^)- Then, the effective magnetic field is slowly 
switched off (h is reduced) while the interaction along x is 
turned on by tuning the relative duration of the evolutions 

under the two Hamiltonians H Isin and H Isin „. If the varia- 
tion of the parameters is sufficiently slow, then the system re- 
mains in the instantaneous ground state. And eventually it ap- 
proaches the ground state of the final modified target Hamilto- 
nian H f ina i - -ti/2 Xj,k Jj,k (°~ Z j cr l + °~} cr t ) wnere tne effective 
magnetic field is zero and both interactions along x and z are 
present. This Hamiltonian exhibits ground state long distance 
entanglement when the end spins are weakly coupled to the 
bulk [|17|-[22h. 

However, in general the typical harmonic trapping poten- 
tial of linear ion traps induces long range interactions with 
maximum couplings at the end of the chain. Thus, in order 
to simulate ground state LDE, the trapping potential has to be 
carefully engineered and the end spins interactions have to be 
made weak. This can be realized with segmented micro-traps 



as discussed in Sec. lIIIAl 



III. ENGINEERING OF SPIN HAMILTONIANS WITH 
TRAPPED IONS 

In this section we study how to manipulate the coupling 
strengths Jj^ and how other kinds of spin Hamiltonians can be 
designed. To be specific, we present detailed calculations for 
an existing micro-structured ion trap 1128112911 . The principles 
used to obtain the concrete results presented in what follows 
are, of course, applicable to other segmented traps as well. 



A. Tailoring the coupling constants in a segmented trap 

In the following we will discuss how to generate the ax- 
ial trapping potential which results in the desired coupling 
pattern. The theory follows [36), where coupling patterns 
were calculated, when ions are held in a micro-structured trap, 
which means, that each single ion or ion chain is located at the 
bottom of an approximately harmonic potential. The separa- 
tion of potentials is large (on the order 100 fi m) and the cou- 
pling between different sites can become impractically small 
(Hz) for the purpose described in this article. Here, we con- 
sider ions held in common closely separated and anharmonic 
wells, and we tune both the harmonic and anharmonic part to 
obtain the desired coupling pattern. 

The effective potential can be written as 
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where rf is the effective rf potential present at an rf power 
level Prf = Po- <Pi is tne dc potential originating from electrode 
i set to the voltage Uj - Uq. The effect of the rf field on the 
axial potential is neglected in the following discussion. 

The potential </>, of a single segment electrode set to a non 
vanishing voltage is bound to be symmetric with respect to 
the center of that electrode for the approximation of harmonic 
trapping and remains ideally symmetric for periodic boundary 
conditions. This leads also to symmetric ion configurations, 
unless the 'history' of the string (e. g. applying intentionally 
asymmetric potential patterns) was used to split the chain off- 
centre. We are in a regime, where the electrodes are much 
wider compared to inter ion separations, so potential or elec- 
tric field cannot independently be chosen at each ion position. 

With a given voltage pattern applied to the electrodes, and 
an initial guess of ion positions, one can calculate equilibrium 
positions by minimizing the total energy. Note that several lo- 
cal minima are possible, corresponding to different numbers 
of ions in wells. In addition, permutations of ion positions 
yield identical total energies. After the equilibrium ion posi- 
tions have been determined, we calculate the normal modes of 
an ion string, the Zeeman shifts of individual ions and the re- 
sulting coupling constants. For all numerical simulations pre- 
sented in the rest of this paper, we assume an axial magnetic 
field gradient of 50 T/m. 



The desired coupling pattern is symmetric (a bulk chain 
weakly coupled to outer ions) so we consider symmetric po- 
tential configurations and ion positions. The symmetry can 
easily be broken at will, either by applying non-symmetric 
voltage patterns, or by loading different numbers of ions per 
site. We create up to three potential wells, and depending on 
wells up to four basic scenarios can be distinguished. 

For small excursion Ax, of ion i from the equilibrium posi- 
tion, the motion of the chain can be decomposed into normal 
modes, which is equivalent to say that the force Fy on ion j 

depends linearly on the excursion pattern Ax as 
F = MAx 
Ftj has two contributions: 

1 . the sum of changes of all Coulomb forces 

2. the change due to a different position in the external 
potential 

So M can be read such that the column vectors are the 
forces on ion j if ion i is moved by a unit length from its 
equilibrium position. 
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Coupling can thus be interpreted in three steps: 

1 . a flip of one ion i would lead to a position change in 
a magnetic gradient. The size of that change depends 
only on the magnetic moment, the gradient, and the lo- 
cal curvature of the trapping potential 

2. the new position changes the equilibrium positions of 
all other ions. For ion j the size of the change is propor- 
tional to the size of the position change of ion i and the 
local potential curvature at the position of ion j 

3. the position change Axj is translated to a change of the 
level splitting Aojj, This change is proportional to Axj 
and the gradient at the position Xj 

In the following, three potential shapes are discussed and 
values for ion positions, qubit level splittings (for zero offset 
field), normal modes and couplings are given. Note that, for 
symmetric potentials and ion positions, the couplings are also 
symmetric, thus 
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1. Coupling in a common well 

When one or more segments are set to a negative voltage, 
and the others are kept at ground, a single trapping zone is cre- 
ated and the potential is predominantly harmonic. The cou- 
pling pattern shows strong next neighbor coupling and also 



substantial coupling beyond next neighbors. Figure Q] shows 
the potential in the vicinity of the center of symmetry together 
with the equilibrium positions and the resulting coupling pat- 
tern. Numerical values for positions, levels splittings, normal 
modes and couplings are given in Tab. U 




-20 20 

position (pm) 



i 12 3 4 

i ,,(|im) -85.0 -14.9 14.9 85.0 

o>,-/2Jr (MHz) -59.5 -10.4 10.4 59.5 

v,/27r(kHz) 36.1 66.6 85.4 85.5 



J tj (kHz) 



1,2 

0.14 



1,3 1,4 2,3 
0.10 0.00 3.93 



TABLE II: Positions x h qubit splittings cou mechanical normal 
modes v, and couplings /y for an ion chain in three approximately 
harmonic axial trapping potential. 



3. Coupling in a single strongly anharmonic well 



FIG. 1 : Axial potential and equilibrium ion positions for a predomi- 
nantly harmonic single well 



i 12 3 4 

x 0J (urn) -37.7 -12.2 12.2 37.7 

«,/27r(MHz) -26.4 -8.5 8.5 26.4 

v,/27r(kHz) 36.1 61.6 83.3 101.8 



J, j (kHz) 



1,2 

2.64 



1,3 

1.85 



1,4 2,3 
1.32 2.67 



TABLE I: Positions x,, qubit splittings W;, normal modes v, and cou- 
plings Jjj for an ion chain in an approximately harmonic axial trap- 
ping potential. 



2. Coupling in three wells 

Trapping ions in three independent wells allows an intu- 
itive approach to generate LDE: the inner ions are confined 
in a common well and couple strongly. The outer ions sit in 
separate wells and thus show only small coupling to the cen- 
ter 'bulk' string (see fig. [2] left image). In fact, separating 
the outer ions might generate outer wells with a steep confine- 
ment, further reducing the coupling to all other ions. This can 
be seen in the resulting coupling pattern (see fig. [2] right im- 
age). Numerical values for positions, levels splittings, normal 
modes and couplings are given in Tab. [II] 




-100 -50 50 100 

position (pm) 



FIG. 2: Axial potential and equilibrium ion positions for ions con- 
fined in three wells 



Making a well strongly anharmonic, substantially alters the 
normal mode spectrum and allows to generate a pattern suit- 
able for creating LDE. Changing the softness of the outer 
wells allows to choose the ratio of the coupling of outer ions 
to their neighbors with respect to the coupling between the 
two center ions in a wide range (see fig. [3}- Note that this 
potential configuration differs only slightly from the situation 
shown before. This indicates a strong dependence of the cou- 
pling on the applied voltages and puts strict requirements on 
voltage stability and accuracy which have to be taken account 
the the design of the voltage supplies ll37ll . Numerical values 
vor positions, levels splittings, normal modes and couplings 
are given in Tab. [Ill] 




50 50 

position (pm) 



FIG. 3: Axial potential and equilibrium ion positions for ions con- 
fined in a single anharmonic well, whose shape generates a coupling 
pattern suitable for creating LDE. The coupling of the inner ions 
compared to the coupling to the end ions can be scaled by choos- 
ing the anharmonicity. 



i 12 3 4 

x ,,(nm) -79.5 -15.1 15.1 79.5 
w,-/27r(MHz) -55.6 -10.5 10.5 55.6 
v,/2;r(kHz) 36.1 52.9 53.0 65.9 



J, j (kHz) 



1,2 
0.14 



1,3 1,4 2,3 
0.10 0.00 3.93 



TABLE III: Positions x,-, qubit splittings (Oi, normal modes v; and 
couplings Jjj for an ion chain in a strongly anharmonic axial trapping 
potential. 



4. A Suitable Potential to Simulate LDE 

A pattern similar to the previous section is used for the sim- 
ulations in the rest of the article. Compared to Sec. IIIIA3I 
all voltages are scaled up by a factor of four, resulting in ap- 
proximately two times stiffer axial trapping frequencies which 
makes trapping more robust and reduces the populations in 
higher vibrational modes. The couplings on the other hand are 
reduced by a factor of 4, slowing down the speed of the sim- 
ulation and thus requiring longer coherence. The long range 
interactions are small, and the end spins are weakly coupled 
giving appropriate starting conditions to generate LDE. Nu- 
merical values for positions, levels splittings, normal modes 
and couplings are given in Tab. [IV] 



x Q ji (nm) 
C0i/2n (MHz) 
y,/2/r (kHz) 

i, J 
Ju (kHz) 



TABLE IV: Positions x,-, qubit splittings w,-, normal modes v; and 
couplings J/j for an ion chain in three approximately harmonic axial 
trapping potential. 
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(9) 



In each time interval in which Q m + 0, the driving frequency 
is close to resonance to a single spin j m , with a small detuning 
h m , that is v,„ = ojj m — h,„. 

The spin dynamics can be further harnessed by controlling 
the external magnetic field. Here we are interested in situa- 
tions in which the spin-spin couplings Jjk are fixed. In partic- 
ular Jjjf does not change when, for example, the direction of 
the B field is inverted. On the other hand in this case the spin 
resonances are affected. Specifically, the difference between 
the spin resonance frequency 

Aj* = ojj - co k 

changes in sign when the direction of the magnetic field is 
inverted. Hence in each step the spin frequency can be defined 

as 



fin) 



COq + T] m 0)0j, 



where coq is an offset frequency, coqj the frequency relative to 
the offset and r\ m — + 1 accounts for two possible direction of 
the magnetic field. 

The system Hamiltonian is 



Hit) = J]e»W 
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B. Spin dynamics with microwave pulses 

The spin dynamics can be manipulated using external mi- 
crowave fields which can drive selectively a given spin by tun- 
ing the driving frequency to the corresponding resonance. The 
hamiltonian for the interaction between the ion spins and the 
driving field in the rotating wave approximation and neglect- 
ing the ion motion (see Sec. lIVBl for an analysis of the effects 
of the motion) takes the form 



H L (f) = -im(t)^{cr^< 
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In general the amplitude Q, the frequency v, and the phase <p 
of the driving field can be time dependent. In particular we 
consider a sequence of step-like driving pulses: We identify 
a set of time instants t„, with m = 0, 1, 2, ■ • • which define a 
corresponding set of time intervals (t m -i,t m ] during which the 
driving Hamiltonian is constant, (Q(f) = Q„„ v(t) - v m and 
ip(t) = <p m for t € (t m -i,t m J). In certain intervals the driving 
field can also be zero. If we define the square-pulse function 



fm(0 - Q(f - fm-l) - &(t - t m ) 



(7) 



with B(t) = for t < and 8(t) = 1 for t > 0, then Eq. © 
takes the form 



H L (t) = Yj^itWfit) 



(8) 
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where H, = | %" os 'a~. 

The time dependence of the driving Hamiltonian H^'(f) 
can be removed in a reference frame rotating at the driving 
field frequency as discussed in App. [A] The dynamics in the 
two representations are related by a unitary and local transfor- 
mation, thus the corresponding entanglement properties are 
equal. 

The system Hamiltonian in the new representation takes the 
form (see App.lAl 
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Here the sum over the time intervals is divided into two sums 
over the intervals in which the driving field is on and off re- 
spectively, and the label j„, indicates the spin which is driven 



quasi resonantly in each time interval. Moreover we have set 
the driving phase <p m to the value <p m - —<f>. with 

m— 1 
Am) V T (m' + l) (m')l . /n\ 
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when Q m =£ 0. 



(14) 



So, in this representation, the spin j m sees an effective mag- 
netic field along the y-axes. 



C. Stroboscopic engineering of the XX spin dynamics 

In order to engineer the dynamics of an XX quantum 
spin model we consider a sequence of driving pulses made 
of 2N + 2 steps (m - 1, ■ • -IN + 2) and identified by the 
corresponding values of the driven spin j m , the direction of 
the magnetic field n m , the laser intensity O,,, and the duration 
of the time step t m - t m -i. During the all sequence the value 
of the detunings is fixed h m - h. The sequence is made 
as follows: after a free evolution (no driving: m — 1 and 
Qi = 0), of time t\ - to - Afi, each spin is driven sequentially 
from the last to the first (such that for m = 2,- ■ -N + 1, 
{jm =N-m + 2, r\ m = -1, Q m = -Sl jm , t m - t m .\ = Stjjj, 

where the actual values of Q,j m and 8tj m are iden- 
tified below); then after another free evolution (no 
driving: m = N + 2 and Q^+2 = 0), of time 
tN+2 - tff+\ - A?2, the spins are driven sequentially 
from the first to the last (such that for m = N + 3, • • • 2N + 2, 
[jm - m - N - 2, t] m = 1, Q. m = Cl jm , t,„ - t m -\ = 6t jm }). Note 
that in the first train of pulses the values of 77,,, and Q m are the 
opposite of that in the second train of pulses. This is realized 
changing the direction of the magnetic field, that realizes the 
transformation j] m — > -77,,,, and adding a phase of n to the 
driving field, that realizes the transformation Q„, — > -Q. m . 

As discussed in App.lBlthe evolution operator at the final 
time t - t2N+2 - Af 1 + A/2 + Yjj $tj, corresponding to this 
sequence, with the parameters in each time interval which sat- 
isfy the relations 
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?Yj A jjJtj = njjt 



with ttj m e ! 



(15) 



]>}m 



can be approximated, in the limit \Ajj m \ s> \Q, m \ » \h m \ , \Jjjc\, 
as 
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where 



/ N 1 
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0.08 
time (s) 



0.16 



0.02 0.04 

time (s) 



FIG. 4: End-to-end concurrence (thick, blue lines) and fidelity (thin, 
red lines) with the instantaneous ground state for a four ion system. 
The dashed line are obtained integrating the time dependent effec- 
tive Hamiltonian in Eq. ( I22t with the time dependence defined in 
Eq. ( 1231 ) (ho = 2n x 0.99kHz)', The solid lines are obtained using 
the Hamiltonian in Eq. Jill ), following the pulse scheme described 
in Sec. lIII Cl and with the stepwise variation of the parameters h and 
a. The spin-spin coupling constant and the spin resonance frequen- 
cies are identified in Sec. MI A 41 Plots (a) and (b) differ for the ve- 
locity of the adiabatic manipulation: in (a) r = 2n x 10Hz in (b) 
r = 2k x 40Hz- The other parameters are: (a) Afi = 64/is and 



{S u ---6t 4 } = jl.87, 1.95, 1.94, 1.87)^, and (b) A?, 
[Su ■ ■ ■ St 4 ) = (1.57, 1.61, 1.61, 1.57) /is 



SOfis and 



The stroboscopic evolution at times nt, with n e N, given by 
the repeated application of this sequence of pulses is described 
by the operator 



tt jjn l a -iHf>- Ah-iH 1 .^ At,\" 

According to the Trotter formula [..cite] 



,-i(H 1+ H2)t 



lim (> 



^-Olitln^-iHitl 

-1 ,,,-1 



")"- 



(18) 



(19) 



and in the limit Afi, A?2 <K \Jjjc\ >\h\ \ we can approximate 



U„ 



where 



-i(H l '\ +aH ( . x) )«A?i 



Af 2 
Afi 
Afi 

P = -r 



\ I sing I sing f 



(20) 



Afi 



f (1 +a)Afi +2 2j<5f; 



(21) 



This result demonstrates that the stroboscopic evolution de- 
scribed above approximate the evolution, at times nt, of a spin 
system with the effective Hamiltonian 

ff*r=j8K^ + "H«J. (22) 



IV. ADIABATIC PREPARATION AND 
STROBOSCOPIC/PULSED DYNAMICS 

The parameters h, a and /3 can be varied adiabatically in 
order to prepare the ground state of an XX hamiltonian as dis- 
cussed in App. [C]The effective external magnetic field h can 




0.08 
time (s 



0.02 0.04 

time (s) 



FIG. 5: As in Fig.|4]with a dephasing rate (a) y = 2nxOAHz and (b) 
y = In x OAHz. 



on the other and the rate is too large (Fig. |4) (b) then the evo- 
lution is no more adiabatic and the system ends up in a state 
which is not the ground state of the final hamiltonian and the 
fidelity is reduced. Nevertheless also in this case, the end-to- 
end entanglement can be very large at certain times meaning 
that the end spins approach a Bell state. Also in this case the 
bulk efficiently mediates an interaction between the end spins 
which entangle them efficiently. 

In all cases the evolution performed with the effective 
Hamiltonian and that obtained via the sequence of pulses are 
very similar meaning that the protocol is faithful. The effi- 
ciency of the second is only slightly reduced but the overall 
behaviors are consistent. 



be varied by adjusting the detuning between the driving field 
and the ion-spin resonance frequencies. On the other hand a, 
and correspondingly/? are varied by controlling the time A?2- 

The system is initialized in the ferromagnetic state with all 
the spins aligned along the z-axes that is the ground state of 
the Ising Hamiltonian (a = 0). The value of h is initially set to 
some value ho in order to remove the ground state degeneracy. 

The parameters a and h are then slowly varied to realize 
the adiabatic preparation of the LDE. In particular during each 
sequence of pulses, that is described in Sec lHICl the values of 
h and a are kept fixed, while they are varied from sequence to 
sequence in order to realize a stepwise approximation of the 
functions (see App.lCb 



a(t) 

h{t) 



l-e~' 

h e- n . 



(23) 



The efficiency of this stepwise adiabatic protocol is ana- 
lyzed numerically by evaluating the evolution corresponding 
to the Hamiltonian dTTb with the time sequence and the param- 
eters discussed in Sec. MI CI and the corresponding stepwise 
variation of h and a. The results are shown in Figs. [4] and [5] 
(solid lines) for different values of the rate of variation r. They 
are compared with that obtained by the numerical integration 
of the Shrodinger equation with the effective time dependent 
Hamiltonian|22]where the time dependent parameters a and h 
are defined in Eq. ( |23l (dashed lines). 

In both cases the protocol is characterized in terms of the 
fidelity between the resulting state and the expected instan- 
taneous ground state of the effective Hamiltonian d22b (red, 
thin curves), and in terms of the end-to-end concurrence (blue, 
thick curves). The fidelity indicates how far is the resulting 
state from the expected one: Fidelity equal to one corresponds 
to perfect adiabatic following; while equal fidelity for both 
the standard adiabatic evolution (dashed lines) and the step- 
wise adiabatic evolution (solid lines) means that the protocol 
realizes a perfect simulation of the effective Hamiltonian. On 
the other end, the concurrence measure the entanglement be- 
tween the end spins, and concurrence equal to one indicates a 
maximally entangled Bell state. 

When r is sufficiently small (Fig. |4| (a) the ground state 
preparation is quite good: The fidelity is close to one and as 
expected the ground state exhibits large entanglement. When 



A. Effect of spin dephasing 

Fig.|5]is evaluated including the dephasing of the spins. It is 
obtained by solving a master equation for the spins dynamics 
of the form 



p = -i[H(t),p]+£ D p 



(24) 



where H(f) corresponds to the Hamiltonian dTTT > for the solid 
lines and to the effective Hamiltonian ( l22l for the dashed lines. 
Moreover _£d accounts for the spins dephasing at rate y and 
takes the form 



£ Dp = lTii a *jp aX j-p)- 



(25) 



As expected, the dephasing reduces the efficiency of the 
scheme. In Fig.[5](b) the entanglement of the end spins is still 
relatively high with a dephasing time of ~ 0.4s. The scheme 
is expected to be resistant to stronger dephasing (smaller de- 
phasing time) when the coupling constants Jj^ are larger. This 
could be obtained with a stronger gradient of the external 
magnetic field. 



B. Mechanical effects 

So far we have neglected the motion of the ions. Internal 
electronic dynamics and motion can be coupled by an electro- 
magnetic filed. In particular when the ions are in a magnetic 
gradient also long wavelength radiation, as microwaves, can 
have a significant mechanical effect allowing for example for 
sideband cooling [30, 31, 39]. In the following we justify our 
treatment in which we neglect the atomic motion. 

In a magnetic gradient the coupling between an ion j and 
a mechanical normal mode k is scaled by the effective Lamb- 
Dicke (LD) parameters 113111 



rij,k 



h n B g dB 



2m LOk h o)k dx 



' U 



(26) 



X-Xqj 



where TJH^ is the frequency of the normal modes, and S is 
the matrix that diagonalize the Hessian matrix A (see SeclIB 
of the potential energy function that confine the ions, that is 



(S T AS) . = 5j,k majj. These parameters are typically small 

and allow for a systematic expansion of the corresponding dy- 
namics in power of 77^. Including the lowest order mechan- 
ical effects the Hamiltonian for the interaction between the 
ions and the driving field (see also Eq. ©) takes the form 



H L (t) = -ino.it) 

x ZW l+ Tj i] j*{ a i- ak ) 



(27) 
■4p»r*rffl] _ ft 



where h.c. stands for the hermitian conjugate, and a k , ak 
are the creation and annihilation operators for the vibrational 
mode k. This Hamiltonian accounts for sideband transitions 
at frequencies Uj + u>k. In order for our analysis to be justi- 
fied the corresponding excitation probability have to be very 
small. It can be estimated as 



O(0 J]j, k 



V(t) - U)j±0Jk 



(n k + 1) <§C 1 



(28) 



where iik indicates the average number of motional excitations 
in the mode k. For the parameters used in Figs. |4]and|5] the 



(max) 



largest LD parameter has a value of 77 . , ~ 0. 14 and the low- 
est vibrational frequency is To 



(inin) 
k 



71 kHz. On the other 



hand Q. - n/(48tj) ~ 70-80 kHz and v(f) - coj is in the 
interval [0, ho], with ho ~ 0.99 KHz, for the spins which are 
driven close to resonance, and in the range [0.07, 0.126] GHz 
for the other spins. The corresponding largest sideband tran- 
sition probability is ~ 0.027 (n^ + 1), which demonstrate the 
validity of our results also for Doppler cooled trapped ions, 
without additional sub Doppler cooling to the ground state of 
the axial potential. Increasing the gradient of the B field, the 
coupling strengths increases allowing for a faster preparation; 
However the system approaches the regime in which the me- 
chanical effects are relevant. In fact, stronger gradient of the 
magnetic field correspond to larger LD parameters. In turns 
it imply that the LD expansion up to the first sidebands is no 
more valid and higher order sideband transition may become 
relevant hence invalidating our treatment. 



C. Experimental Feasibility 

Here, we first consider how a fluctuating coupling constant 
(due to, for instance, fluctuating electrode voltages) would in- 
fluence the generation of LDE. Then we briefly outline typical 
experimental parameters, as are used in the experimental set- 
up described in 112911 . 




FIG. 6: End-to-end concurrence for random couplings: for each 
value of the interval of variation A£ of the random variables f jJr the 
plot shows the concurrence corresponding to 500 different realiza- 
tions (the 500 points form a vertical segment for each value of A£). 
The horizontal line indicates the value 0.92, that is the value of the 
concurrence for the coupling matrix defined in Sec. IIII A 41 



Let us now define the random couplings Jj^ - Jjk ( 1 + %jjA 
with £,£ uniformly distributed random variables in the range 
[-A£, A£]. The resulting concurrence are reported in Fig. |6] 
These results shows that the end-to-end entanglement is suffi- 
ciently stable under variations of the spin-spin couplings. 

In 112911 the operation of a segmented trap with a built-in 
switchable gradient based on micro-structured solenoid is re- 
ported. Different isotopes of Ytterbium with or without hy- 
perfine structure can be trapped. For the experiment proposed 
here, we use 171 Ybwith a nuclear spin of / = 1/2 yield- 
ing two hyperfine levels with F — 0, F — 1 in the elec- 
tronic ground state 113 3J 14311 . Different qubit implementations 
are possible and magnetic sensitive states can be included 
to allow for magnetic gradient induced coupling (MAGIC) 
II30h32I1 as required for the experiments discussed here or 
avoided to yield a quantum memory with a long coherence 
time. The qubit is manipulated using microwave fields (ap- 
plied through a conventional wave guide) and Rabi frequen- 
cies near Q a; 100 kHz on the cr transition and a bare coher- 
ence time of the magnetic field sensitive hyperfine qubit on 
the order of 5 ms have been observed. Applying spin echo 
techniques [40], dynamical decoupling [41] or dressed states 
B42I1 . we expect to be able to observe a coherent time evo- 
lution on a second timescale. Splitting and merging of ion 
strings (which involves the generation of anharmonic trapping 
potentials), as well as shuttling have been demonstrated. Sta- 
ble trapping could be observed down to axial trap frequencies 
around v\ — 2n ■ 40 kHz. 



FIG. 7: Chain of 33 Yb ions in our segmented trap. 



1. Stability estimates 

Let us consider a spin model with Hamiltonian H = 

2/i Jj,k l "/ ^ + crZa k)- ^ tne cou pli n g matrix J — J where J 
is the one defined in Sec. IIII A"4l then the ground state end-to- 
end concurrence is 0.98. 



V. CONCLUSION AND OUTLOOK 

In this theoretical article we have introduced and inves- 
tigated schemes for the quantum simulation of LDE with 
trapped ions. In particular, we have shown how to tailor the 



trapping potential in order to engineer a specific spin-spin 
coupling pattern in one-dimensional lattices, and we have de- 
signed a sequence of microwave pulses able to engineer ef- 
fective spin- 1/2 Hamiltonians of the gapless and isotropic XX 
type. Finally, we have analyzed the efficiency of the adia- 
batic quantum preparation of the ground state of an effective 
Hamiltonian which exhibits LDE, demonstrating its feasibility 
within the limits of current ion trap technology. In the course 
of the investigation, we have introduced and combined for the 
first time trap shaping, adiabatic preparation, and Trotteriza- 
tion of the interactions, three key elements that have never be- 
ing combined together and that we have proven necessary for 
the realization, so far not yet attained, of truly and highly non- 
classical features of complex models of interacting quantum 
many-body systems. 

We expect to be able to carry out the experiment in a 
medium-term time scale as detailed in IIVCI To improve the 
level of control over the anharmonicity of the axial trapping 
potential, it might be necessary to use smaller axial trapping 
segments, possibly in a surface trap. Larger gradients would 
boost the coupling and allow for stiffer axial confinement, 
making the scheme more robust against thermal excitation and 
external stray fields. 
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Appendix A: The dynamics in a new reference frame 

Let us consider the model described by the Hamiltonian in 
Eq. ( TTOl ). and study the dynamics in a new reference frame 
defined by the unitary transformation 



U (t) = ^U^^'^'V 1 ^ 



'(t,„-l-t,„-2) ... e -i<'(fl-fo) 



+0(t - 1) 



with 



H, 



(«!) 



2 ^-J J 



"V. 

J 



(Al) 



(A2) 



where the actual values of v are identified below. Since the 
unitary transformation is local, the entanglement properties in 
the new representation are the same as that in the original one. 
If \iff(t)) is the state in the original representation, then the 
dynamics of the transformed state \i[/(t)) - Ul(f)\if/(f)) is ruled 



-(B)) 



by the Hamiltonian 

H(t) = ul(t)H(t)U (t) - J] e m {i)Ht 

m 

= ^e m (0[^ m) +// zz + 4 ffl) ] (A3) 

m 

with 

H (,n) = gO") _ H 0n) = 1 ^{Jf ~ V$»V* 

J 

Hf' } = -iO m ^ L+^-^'X'^A - h.c\ (A4) 
i 
where 

m— 1 

Am) V" 1 f (m' + l) On')']. / a c\ 

m' = l 

The last Hamiltonian is obtained exploiting the relation 
e 1 * Arte -1 * ■*' = <r|e 21 ^'. The values of v are chosen so 
that the new Hamiltonian is time independent: 

y (m) _ Jm) _ ^ when Dm _ q 

y (m) _ Jjn) _ ^ when Dm + q_ (A6) 

Moreover the phase is fixed to the value 



<Pm - ~<Pj m 



(Al) 



Am) 



so that the phase of the driving pulse, (jr. + ip m , is zero for 
the spin j m which is the one that is close to resonance with 
the driving field in each time interval m. Thereby we obtain 
Eq. CD). 

The corresponding evolution operator decomposes into the 
product of evolution operators each one describing the evolu- 
tion in a single time interval [t m -i,t m ]: 

U(t) = J] e m (0e- iH, "' ,( '-'"'-' ) e- iH, '""' ( '"-'-'"- 2) • ..^ { "("-'«\ 



where 



H im) = H (m) + ^ + H (m)_ 



Appendix B: The sequence of driving pulses 



(A8) 
(A9) 



We are interested in the limit in which \u>j - w Jm | » |Q m | » 

\hm\ > \Jj,k\, for j + jm- Hence we can approximate the Hamil- 
tonian ( TTlT i by retaining only the resonant terms as 



H(t) - Y, *M[h$> +h «] 

m : n,„=0 

+ J] e m (t)[H^ + Hf"] 

m : n,„#0 



H 



(m) 



hQ m o y jm 



(Bl) 
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The evolution operator corresponding to the sequence of 
pulses described in Sec. |HI CI can be written, using the ap- 
proximate Hamiltonian (IB II) . in the form 



■j(+)a 



i// 1 , .«. . . e -iHT'<»i e -i#S„ s A'2 e -iH\->6h . . . e -iH' N >Sl N 

-in, . za/1 






(B2) 



where 



7=1 j* 

QlO^ + - V A/ j oi 



j*jm 



and the total time of the sequence is 



t - At\ + Ah + 



!>• 



The operator Uj can be rewritten as 



where 



with 



Uf = T N ---Tie '■""* -T[...T N t '™« 



r /n? = e ~ io "»' T ie~ in '»<' r M ,5 '-"» 



£>7,„ = ^ X Aju.fy 



7>7». 



Now we use the relation 
S (T, A) = e - lA ^ e - ir ^o^e ir ^e iA ^ 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



= cos (2T) cxj + sin (2T) [cos (2A) cr) + sin (2A) cr^.] , 

which reduces to S (F, A) = crj when r = \+nn and A = n'n, 
or when r = — ? + wr and A = ±| + n'7r with n, n' € Z. Thus 



setting, for example, 



%A, 



£)/,„ = W/^TT 



(B9) 



with rtj m e Z, V/ m , then 



„ -i//< v > Af 2 -i/f' :) . Af, 



(BIO) 



where 



, N . 

7=1 j* 



0.0 0.5 1.0 



Gap 



-(kHz) 




0.5 
h/k, 0.0 



^•ci 




(t) 



(h) 



h(t)/h 



1/r 5/r 10/r 



0.5 



(c) 



0^— - 



0.08 
time (s) 



0.16 



FIG. 8: (a) Gap between ground a first excited state of H^ in Eq. <22t 
with four ions, in the space of parameters {h, a}. The red line in- 
dicates the gap corresponding to the adiabatic variation of a and 
h. (b) time evolution of the parameters a{t) and h{t). (c) End-to- 
end concurrence (thick, blue line) and fidelity with the instantaneous 
ground state (thin, red line), obtained integrating the time dependent 
Shrodinger equation with the effective Hamiltonian in Eq. J22I ) with 
four ions, and with h = 2n x 0.99kHz and r = 2n x lOHz- The 
spin-spin couplings are defined in Sec. lIII A"4l Although not relevant 
for the present result, in order to be consistent with the results of 
Sec. lIVI we have set the parameter /3 to the values defined in Eq. i2lt 
(a different value of /3 corresponds to a rescaling of the energy and 
correspondingly of the duration of the protocol). 



Appendix C: Adiabatic preparation of the Ground state 

A system initially in an eigenstate |i///(0)) of its Hamilto- 
nian, follows the instantaneous eigenstate \iffj(t)), which de- 
rive from the initial state by continuity, when the correspond- 
ing Hamiltonian is deformed adiabatically [44]. Condition for 
the adiabatic evolution is that during the evolution the proba- 
bility for the transition form the eigenstate \i/fj(f)) to a different 
one |tAt(0) (V&) is negligible, this can be estimated as [44] 



I 



mt)\dH(t)ldt\^j(t)) 



[E k (t) - Ej(t)] 



« 1. 



(CI) 



It means that larger is the difference in energy between the 
eigenstate state |i///0) and all the other, more easily the adia- 
batic condition can be satisfied. 

In particular if initially the system is prepared in the ground 
state then it will remain in the instantaneous ground state un- 
der a slow variation of some Hamiltonian parameters. This 
idea can be applied to prepare the ground state of compli- 
cated Hamiltonians: One can first prepare the ground state 
of a sufficiently simple one which is easy to prepare. Then the 
Hamiltonian is adiabatically changed until approaching the fi- 
nal target Hamiltonian. Correspondingly the system will end 
up in the ground state of the final Hamiltonian. 

In our case according to the result of Sec. |IH CI we are able 
to generate the dynamics corresponding to the Hamiltonian 
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(see Eq. d22i i) 



i 

-Pit) J] Jj, k [0)0% + a(t) cry k ] . (C2) 
;■* 

where /3(f) is function of a(t) as specified in Eq. (TSTT i. 

We want to prepare the ground state of H X z = H e ff(l,0) = 
—ft Yijjk Jj,k f ^ ^ + cr / cr |)' Hence we can first prepare the 
ground state of a ferromagnetic Ising Hamiltonian i/i S i ng = 
H ef f (0, /zo) = P (ho Zj Vj ~ Zj,k J i,kO- Z j<T z j) which simply cor- 
responds to the ferromagnetic state in which all the spins are 
polarized along z. Then the ground state of Hxz is obtained 
by the adiabatic variation of the parameters hjho : 1 — > and 
a: 0->l. 

An example of adiabatic preparation of the ground state of 
the Hamiltonian H e g ( 1 , 0) is shown in Fig [8] The parameters 
a and h are varied according to (see the curves Fig. |8](b)) 



h(t) = fc e" rt 
a(0 = 1 - e" 71 



Initially the parameters can be varied rapidly because the cor- 
responding gap between ground and first excited state is rela- 
tively large as depicted in Fig. [8] (a). As H^ approaches the 
target Hamiltonian the gap reduces and correspondingly the 
variation have to slow down. The curves in Fig.|8](c) are ob- 
tained by numerical integration of the Schrodinger Equation 
with the time dependent Hamiltonian dC2l ). and are equal to 
the dashed curves in Fig.|4](a). The red, thin line in Fig.|8](c) is 
the fidelity between the state obtained with the adiabatic evo- 
lution and the instantaneous ground state. This curve is very 
close to 1 at all times indicating that the system actually fol- 
lows the adiabatic ground state. The spin-spin couplings that 
are used in these calculation are that discussed in Sec. lIII A"4l 
The Hamiltonian H e g(l,0) with these coupling strengths ex- 
hibits long range entanglement, that is strong entanglement 
between first and last spin. This feature is described by the 
blue tick curve in Fig.|8](c), that displays the entanglement, as 
measured by the concurrence between first and last spins. As 
expected, at large time the end spins are strongly entangled. 
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